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Abstract 

When averages of different experimental determinations of the same quantity are computed, 
each with statistical and systematic error components, then frequently the statistical and sys¬ 
tematic components of the combined error are quoted explicitly. These are important pieces of 
information since statistical errors scale differently and often more favorably with the sample 
size than most systematical or theoretical errors. In this communication we describe a trans¬ 
parent procedure by which the statistical and systematic error components of the combination 
uncertainty can be obtained. We develop a general method and derive a general formula for 
the case of Gaussian errors with or without correlations. The method can easily be applied to 
other error distributions, as well. For the case of two measurements, we also define disparity 
and misalignment angles, and discuss their relation to the combination weight factors. 


1 Introduction 

Error propagation as well as the averaging of results of individual measurements — at least in 
the context of strictly Gaussian errors, possibly with statistical or systematic correlations — are 
straightforward, they are covered in many textbooks^ and there seems to be no open issue, because 
all that is required is multi-variate analysis applied to normal distributions. It is the more surprising, 
that — to the best of the authors knowledge — no explicit analytical expression is available that 
serves to compute for a given set of measurements of some quantity with individual, generally 
correlated, errors of statistical and systematic nature, the statistical and systematic components of 
the uncertainty of the average. 

That is to say, while in the Gaussian context it is clear how to obtain the average including its 
uncertainty, and that the total error ought to be the quadratic sum of the statistical and systematic 
(and perhaps other such as theoretical) error components, formulae of these individual components 
and their derivation have not received much attention. But the proper disclosure of the statistical 
(random) error component compared to the systematic uncertainty can be of importance. For 
example, in the context of the design of future experimental facilities it is crucial to know how much 
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1 For a readable and practical treatise from the Bayesian point of view see Ref. [T]. This reference strikes a welcome 
balance between theoretical background and applied data analysis offering more than recipes. For an upcoming 
textbook on data analysis in Particle Physics, see Ref. [2]. 
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precision (say, over the world average) can be gained by simply generating larger data samples, in 
contrast to possible technological or scientific breakthroughs. The systematic error is, of course, 
more troublesome as it cannot be reducecfl as straightforwardly by increasing the sample size N. 

In the next Section, we review the standard procedure to average a number of experimental 
determinations of some observable quantity. We also mention an approximate method to obtain the 
statistical and systematic error of an average of similar experiments where the ratio of the systematic 
to statistical components are comparable, or where the statistical error is dominant. 

Then we turn to the main point, the exact determination of the individual error components of 
an average. We show that in the absence of correlations the square of the statistical error or any 
other type of uncertainty is weighted by the fourth power of the total errors. 

Then we turn to correlations, starting with the simplest case of two measurements for which 
we introduce the concepts of disparity and misalignment angles. Finally, we present exact relations 
for the case of more than two measurements, and address some problems that arise when new 
measurements are added to an existing average iteratively. 


2 Simplified Procedures 

Suppose one is given a set of measurements of some quantity v. with central values Vi, statistical 
(random) errors r,; and total systematic errors Si. For simplicity, we are going to assume that the r* 
and Si are Gaussian distributed (the generalization to other error distributions is straightforward) 
in which case the total errors of the individual measurements are given by 


U = 


rf + sf 


(1) 


If we furthermore temporarily assume that the measurements are uncorrelated, then the central 
value v of their combination is given by the precision weighted average, 


v = 


E i Vjtj 

E i E 


-2 


with total error 



Similarly, the statistical component r of t can often be approximated by 


( 2 ) 

(3) 


The systematic component s of t is then obtained from 

s = sjt 2 - f 2 . (5) 

2 Simply collecting more data often helps to reduce even the systematic error component, because some error 
sources that are nominally classified as systematic may be tracable to represent themselves statistical measurements. 
Moreover, with growing statistics one may restrict oneself to cleaner data by imposing stronger selection criteria 
(cuts). Nevertheless, the N ~ 1 / 2 scaling of the statistical component is usually out of reach. 


(4) 
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experiment 

Vi 

Ti 

Si 

U 

Ui 

ALEPH 

0.1451 

0.0052 

0.0029 

0.0060 

0.0057 

DELPHI 

0.1359 

0.0079 

0.0055 

0.0096 

0.0095 

L3 

0.1476 

0.0088 

0.0062 

0.0108 

0.0106 

OPAL 

0.1456 

0.0076 

0.0057 

0.0095 

0.0094 

combination 

0.1439 

0.0035 

0.0026 

0.0043 

0.0040 


Table 1: Central value, Vi, random error, r,;, systematic uncertainty, Si, total error, £,;, and un¬ 
correlated error component, Ui, for measurements of the quantity A T by each of the four LEP 
collaborations, as well as the corresponding numbers for the combined result. 


For example, two measurements with 


ri = S 2 = 30, V2 = si = 40, £i = £2 = 50, (6) 


would result in 

r = 24, s = V674 w26« r. (7) 

Notice, that while the individual errors in Eq. ([6]) are symmetric under the simultaneous exchange 
of the statistical and systematic errors (we recall that all r* and .s, are assumed Gaussian) and the 
labels of the two measurements, the result 0 does not exhibit the corresponding symmetry which 
would imply f = s exactly. The exact procedure introduced in the next section would indeed yield 
f = s in this example. 

Now consider the case where one of the systematic errors, say si, is larger and eventually 
si —> 00 . Then the weight of the first measurement approaches zero, and t — > £ 2 , as expected. 
However, one would also expect that f — > and s —> S 2 , while instead r < remains constant 
and s —> \/1924 44 > S 2 - Thus, one would face the unreasonable result that averaging some 

measurement with an irrelevant constraint (with infinite uncertainty) will decrease (increase) the 
statistical (systematic) error component, leaving only the total error invariant. In other words, if 
in a set of measurements there is one with negligible statistical error, then the average would also 
have vanishing statistical error, regardless of how unimportant that one measurement is compared 
to the others. Clearly, Eq. 0 is then unsuitable even as an approximation. 

One can easily extend these consideration to the case where the individual measurements have a 
common contribution c entering the systematic error. The precision weighted average 0 and total 
error 0 are then to be replaced by 


v = 


E i 


t = 


1 


Ei u i 2 V E* u i 

where the uncorrelated error components are given by 

Ui = Jt ? - c 2 , u = \/t 2 - c 2 = 


E; u i 


-2 


( 8 ) 

( 9 ) 


and where tj > 0 requires (? > —u\ for all i. 
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channel 

Vi 

Vi 

Si 

U 

Ui 

CC electron 
CF electron 

0.2302 

0.2312 

0.0009 

0.0007 

0.0013 

0.0013 

0.0016 

0.0015 

0.0012 

0.0011 

all electron 

muon 

0.2308 

0.2307 

0.0006 

0.0009 

0.0012 

0.0012 

0.0013 

0.0015 

0.0008 

0.0011 

all lepton 

0.2307 

0.0005 

0.0011 

0.0012 

0.0007 


Table 2: Same as Table [T| but for the weak mixing angle determinations by ATLAS. 


The general case of correlated errors will be dealt with later, but we note that the case of two 
measurements with Pearson’s correlation coefficient p can always be brought to this form with c 2 
given by 

c 2 = f 2 - u 2 = ptitz = p \Ju\ + c 2 + c 2 . (10) 

A proper (normalizable) probability distribution requires \p\ < 1, so that from Eq. (HOD . 


c 2 > -- 


( 11 ) 


guaranteeing that t is real. On the other hand, iti or 1 x 2 , as well as u , in Eq. © may become 
imaginary provided that 


P 



or 



( 12 ) 


in which case the first or second measurement, respectively, contributes with negative weight, and 
v lays no longer between vi and V 2 - In this situation, one rather (but equivalently) regards the 
measurement with a negative weight as a measurement of some nuisance parameter related to c. 
Replacing the inequalities m by equalities, gives rise to an infinite weight (one of the iq = 0) as 
well as u = 0 and t = c. 

As a concrete example, each of the four experimental collaborations at the Z boson factory 
LEP 1 0 [H 0 Ej has measured some quantity A T (related to the polarization of final-state r leptons 
produced in Z decays) with the results shown in the Table below. A number of uncertainties affected 
the four measurements in a similar way, leading to a relatively weak correlation matrix [2] which, 
while not quite corresponding to the form ©, ©, can be well approximated by it when using the 
average of the square root of the off-diagonal entries of the covariance matrix c ss 0.0016. 

The values in the last line are v, r, s, t and u as calculated from Eqs. (©, ©, © and ©. v, r 
and s agree with Table 4.3 and t agrees with Eq. (4.9) of the LEP combination in Ref. [7]. 

Table [2] shows the more recent example of the determination of the weak mixing angle [8] which 
is based on purely central (CC) electron events, events with a forward electron (CF), as well as muon 
pairs. Here the average of the off-diagonal entries of the covariance matrix amounts to c ~ 0.0010. 
This is an example where the dominant uncertainty is from common systematics, namely from the 
imperfectly known parton distribution functions affecting the three channels in very similar ways. 

We will return to these examples after deriving exact alternatives to formula ©). 
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3 Derivatively Weighted Errors 


Our starting point is the basic property of a statistical error to scale as N - 1 / 2 with the sample size. 
To implement this, we rewrite Eq. m as 


U = V e2r i + S 1 

Thus, the statistical error satisfies the relation, 


6=1 


,, dU 

n = V (> * 

In the absence of correlations we can use Eq. ©, and demand that analogously, 


_ 9 - dt 

r =t Je 


, dt 




p 4 


(13) 


(14) 


(15) 


Notice, that Eq. © can be recovered from Eq. (ITSl) upon substituting U —> Ti and t —> r. Eq. m 
means that the relative statistical error of the combination, x, is given by the precision weighted 
average 


x 2 = 


2,-2 


E,t - 2 


(16) 


where 


Xi = — 
I'i 


r 

x = - 
t 


Furthermore, giving the systematic components a similar treatment, we find 

T 4 

M, 


s_2 = E s 'i - 


(17) 


(18) 


so that the expected symmetry between the two types of uncertainty becomes manifest, and more¬ 
over, Eq. © now follows directly from Eqs. (usd and m, rather than being enforced. The central 
result is that for uncorrelated errors, the squares of the statistical and systematic components (or 
those of any other type) of an average are the corresponding individual squares weighted by the 
inverse of the fourth power of the individual total errors, or equivalently, weighted by the square of 
the individual precisions t~ 2 . 

Returning to the case where the only source of correlation is a common contribution c ^ 0 
equally affecting all measurements, we find from Eq. (151) . 


f 2 




where 



Applied to the case of A T measurements we now find 


-2 = Ez Vi «i 


2„,—2 


Ei 


-2 



(19) 

( 20 ) 


r = 0.0035, s=VP-r 2 = 0.0026, (21) 

which agree not exactly, but within round-off precision with the approximate numbers in Tabic [T| 
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4 Bivariate Error Distributions 


As a preparation for the most general case of N measurements with arbitrary correlation coefficients, 
we first discuss in some detail the case N = 2. Recall that the covariance matrix in this case reads 


T = f I 1 

K pt\l2 


phh\ _ (tj c 2 

C 2 to 


2 = 


The precision weighted average is given by the expression, 

V\t 2 + V 2 t\ — (Ul + V2)c 2 V1+UJV2 


t\ + i 2 — 2c 2 


where 


t{ - c 2 

u> = 


ti — c 2 


c 2 = 


1 + u> 
t\ — wf 2 

1 — OJ 


obtained by minimizing the likelihood following a bivariate Gaussian distribution, 

L oc e~ x 1 2 , 

where 

X = v 1 v , v = 

The one standard deviation total error t is defined by 

A\' 2 = X 2 (v + t)~ X 2 (v) = 1, 

which results in 


Vi — v 

V2~V 


t = 


/ 1_p2 

/ t\p 2 - C 4 

jtl-uHl 

/fr 2 + * 2 “ 2 -2ptrV ' 

\] t\ + t 2 - 2c 2 \ 

1 1-CJ 2 


or conversely, 


\t\ - P 
t\-P 




( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 


Eg. (1291) is useful in practice if one needs to recover the correlation between a pair of measurement 
uncertainties and their combination error. 

We now turn to the generalization of Eq. m in the presence of a systematic correlation. When 
applying our method of derivatively weighted errors to Eq. (l28l) it is important to keep c 2 = ptit 2 
fixed (this would be different in the presence of a statistical correlation). Doing this, we obtain 


r = 


\frl{t 2 2 -c 2 ) 2 +rl{t\ -c 2 ) 2 v^i + w2r : 


t\ + t\ — 2c 2 
For the systematic component we find 


1 + OJ 


s = 


1 /s 2 + 2wc 2 + w 2 s| 

1 + OJ 


(30) 


(31) 
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and we also note that 


u 2 = - -(32) 

1 — or 

More generally, one can compute the error contribution q of any individual source of uncertainty 
q to the total error as 

J Qi + 2 uc 1 + u 2 ql 

-T+U-’ (33) 

where c 2 is the contribution of q to c 2 with the constraint 

£c 2 = c 2 . ( 34 ) 

Q 

If the two uncertainties q,; are fully correlated or anti-correlated between the two measurements, 
then 

2 , - <?i±wg 2 

c q = ±qiq 2 , q = ——- , (35) 

where the minus sign corresponds to anti-correlation. 

The formalism is now general enough to allow statistical correlations, as well. As we will illus¬ 
trate later, knowing all the q is particularly useful if one wishes to successively include additional 
measurements to a combination one by one — rather than having to deal with a multi-dimensional 
covariance matrix. This situation frequently arises in historical contexts when new measurements 
add information to a set of older ones, rather than superseding them. But there is a problematic 
issue with this, which apparently is not widely appreciated. 

5 Disparity and Misalignment Angles 

Continuing with the case of two measurements, we can relate p to the rotation angle necessary to 
diagonalize the matrix T. If we define an angle 0 quantifying the disparity of the total errors of two 
measurements through 

tan f— = ——— , (36) 

2 ti + 1 2 

then the diagonal from of T is RTR T with 


COS \ 


tana = pcot /3, 


where 


-< a < 

2 ~ ~ 




The angle a may be interpreted as a measure of the misalignment of the two measurements 
with respect to the primary observable of interest v. Uncorrelated measurements of v are aligned 
(p = a = 0), while the case \p\ | tan/3| is reflective of a high degree of misalignment. Indeed, in 

the extreme case where 0 = 0 (|a| = 90°) two correlated measurements (p ^ 0) of the same quantity 
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v are equivalent to two uncorrelated measurements, only one of which having any sensitivity to v at 
all. To reach the decorrelated configuration involves subtle cancellations between correlations and 
anti-correlations of the statistical and systematic error components of the original measurements. 

We can now express the weight factor u in terms of the disparity and misalignment angles P 
and a, 

1 + sin /3(1 — tana) 1 + sin/3 — pcos/3 

(jj — - — ——--—- ( 4(j j 

1 — sin (3(1 + tana) 1 — sin 0 — p cos P 
In the case p = a = 0 this reduces to 


oj = tan 


2 




(41) 


and Eq. (E^l) now reads 


v = 


Vl + V2 
2 


— sin /? 



(42) 


One can write equations of the form m and SI for p ^ 0, as well, with shifted angles P related 
to P by 

esc P = esc P — tan a . (43) 

However, this ceases to work out in the presence of a negative weight (w < 0), in which case one 
would need to replace the trigonometric by the hyperbolic functions. 


6 Multivariate Error Distributions 


To treat cases of more than two measurements with generic correlations, one can choose one of 
two strategies. Either one effectively reduces the procedure to cases of just two measurements (in 
general at the price of some precision loss) by iteratively including additional measurements, or one 
deals with a multi-dimensional covariance matrix. 

We first discuss the latter approach, starting with the trivariate case where 


_( 

T = P3tlt2 


The average can be written as 


with 


P3tlt2 

p2tit 3 \ (i\ 

c 2 

c 3 


t 2 

l 2 

pit 2 t 3 = \4 

t 2 

cl) 

Plt2t3 

4 / V c 2 

r 2 

C 1 

til 


+ L 0 2 v 2 + U 3 V 3 



UJ\+U)2+ W3 



l)(*3- 

C 2 )-( C 2 - C 2)( C 2 

— c 

g), 


u2 = (tl-4)(4-cl)-(4-cl)(4-4), 


The total error is given by 


t = 


det T 


Wl + W 2 + W 3 


(44) 


(45) 

(46) 

(47) 

(48) 

(49) 
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and for its statistical and systematic components we find (in the absence of statistical correlations), 


\/w\r\ + w\r\ + wfrf 
UJ\+U)2+ W3 

respectively. The generalization of Eq. ( 1551 ) is now also straightforward. E.g., in the case of 100% 
correlation between the three measurements we have, 


s = 


\jY,i “hi + “iUjTij 
U>1 + 0J2 + W3 


(50) 


uiqi +^ 2 ff2 +u 3 q 3 
U>1 + W2 + W 3 


(51) 


Analogous expressions hold for cases with N > 3 measurements. For example, the covariance 
matrix for the case of N = 4 reads 


(t\ 

r 2 

c 12 

r 2 

c 13 

cf 4 \ 

C 12 

t 2 

l 2 

C 2 

c 23 

c 24 

c 13 

c 23 

*3 

c 34 

w 4 

C24 

C 34 

t\) 


All that remains to be computed are the weight factors cc,;. We found a convenient expression for 
them, e.g., 


t 2 - r 2 
c 2 c 12 

r 2 — r 2 
c 23 c 12 

r 2 r 2 

c 24 c 12 


r 2 -r 2 
c 23 c 13 

t 2 - r 2 
l 3 c 13 

r 2 - r 2 
c 34 c 13 

(53) 

r 2 r 2 

c 24 c 14 

r 2 — r 2 
C 34 c 14 

f 2 r 2 

£4 c 14 



Thus, the u>i can be obtained by computing the determinant of a matrix which is constructed by 
subtracting the i th column from each of the other columns (or the i th row from each of the other 
rows) and then removing the i th row and column. The reader is now equipped to handle cases of 
any N exactly. 

The alternative strategy to compute averages is to add more measurements iteratively. We illus¬ 
trate this using the example of the ATLAS results on the weak mixing angle (see Table [2j. Besides 
the statistical error (there were no statistical correlations) there were seven sources of systematics, 
six of which correlated between at least two channels. The breakdown of these uncertainties as 
quoted by the ATLAS Collaboration [Sj are shown in Table [5] 

Good agreement is observed with Ref. [8], where the small differences are consistent with pre¬ 
cision loss due to rounding. Indeed, the fact that there are round-off issues can already be seen 
from Ref. [8j, where the quoted total error of the CF electron channel is smaller than the sum (in 
quadrature) of the statistical, PDF, and other systematic errors. A similar issue can be observed 
regarding the quoted combined systematic error which is larger than the sum in quadrature of its 
components. 

However, there are small differences between the results from the exact procedure using Eqs. ((501) 
and m and the iterative strategy using Eq. m and ( 1551 ) . The reason can be traced to the 
asymmetric way in which the error due to higher orders enters the two electron channels. This 
induces subtle dependences of all sources of uncertainties (even those that were initially uncorrelated) 
on the correlated ones. It even affects the uncertainty induced by the finite muon energy resolution, 
which does not enter the electron channels at all. To account for this one can introduce additional 
contributions Ac 2 to the off-diagonal entry of the bivariate covariance matrix of the all-electron 
result and the muon channel. These Ac 2 can be chosen to enforce the exact result, but it is 
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CC 

CF e" 

all e 


Eq. (1331) 

Eq. ([51]) 

Ref. 0 

ACg 

central value 

0.2302 

0.2312 

0.23076 

0.2307 

0.23074 

0.23075 

0.2308 

+46.490 

statistics 

9 

7 

5.6 

9 

4.803 

4.795 

5 

-0.156 

MC statistics 

5 

2 

2.5 

5 

2.376 

2.357 

2 

-0.194 

E e 

4 

6 

3.8 

— 

2.461 

2.490 

3 

+0.308 

A E e 

4 

5 

3.3 

— 

2.144 

2.162 

2 

+0.164 

En 

— 

— 

— 

5 

1.763 

1.764 

2 

+0.010 

PDF 

10 

10 

10 

9 

9.647 

9.647 

9 

-0.011 

higher orders 

3 

1 

1.9 

3 

2.272 

2.255 

2 

-0.161 

other sources 

1 

1 

1 

2 

1.353 

1.353 

2 

+0.001 

total 

15.7 

14.7 

12.9 

15.0 

11.939 

11.938 

11.6 

-0.038 


Table 3: Central values and breakdown of uncertainties (x 10 4 ) of the weak mixing angle determina¬ 
tions by ATLAS. E e and E M refer to the and energy scales, respectively, while A E e denotes 
the electron energy resolution. The last three uncertainties from PDFs, missing higher order correc¬ 
tions, and other sources are taken as fully correlated, whereas the other uncertainties are assumed 
uncorrelated. The fourth column is the average of the two electron channels displayed in the second 
and third columns. The 6 th column adds the muon channel ( 5 th column) to them. The 7 th column 
shows the exact combination result of the three channels. The 8 th column are the corresponding 
numbers quoted by ATLAS. The interpretation of the last column is explained in the text. 


impossible to compute them beforehand. In fact, they depend on the new measurement to be added 
(here the muon channel) and not just the initial measurements (here the two electron channels). 
Moreover, the Ac^ necessary to enforce the correct average central value v differs strongly from 
the Ac^ necessary to enforce the total error t. This observation is a reflection of the fact that the 
combination principle can be violated EH, which we state as the requirement that the combination 
of a number of measurements must not depend on the order in which they added to the average. 

Thus, the iterative procedure generally suffers from a loss of precision. In this example the 
procedure gives nevertheless an excellent approximation because the uncertainty from higher order 
corrections (the origin of the asymmetric uncertainty) is itself very small. But there are cases in 
which the iterative procedure does not provide even a crude approximation and where one should 
use — if possible — the exact method based on the full covariance matrix. Unfortunately, its 
construction is not always possible, e.g ., due to incomplete documentation of past results. Recent 
discussions of related aspects of this conundrum can be found in Refs. nassj. 

7 Summary and Conclusions 

In summary, we have introduced a formalism (derivatively weighted errors) to derive formulas for 
random errors or any error type of uncorrelated Gaussian nature. We introduced what we call 
the disparity and misalignment angles to describe the case of two measurements, and showed their 
relation to the statistical weight factors. For the case of more than two measurements with known 
covariance matrix, we derived some explicit formulas in a form which (as far as we are aware) did 
not appear before. 

It is remarkable, that even in the context of purely Gaussian errors and perfectly known correla- 
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tions there are intractable problems at the most fundamental statistical level. Specifically, they may 
arise even when a number of observations of the same quantity is combined and the error sources 
are recorded and the assumptions regarding their correlations are spelled out carefully. In statistical 
terms, one can conclude that such a combination — despite of all its recorded details — represents 
an insufficient statistics of the available information. The inclusion of further observations of the 
same quantity is then in general ambiguous. 

On the other hand, there is no ambiguity in the absence of correlations or when any correlation 
is common to the set of observations to be combined. The fact that the ambiguities disappear in 
certain limits then reopens the possibility of useful approximations. For example, if an iterative 
procedure has to be chosen, one should first combine measurements where the dominant correlation 
is given approximately by a common contribution. Similarly, the measurements with small or no 
correlation with the other ones, are ideally kept for last. 
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